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ABSTRACT 

We present a deep Chandra X-ray Observatory study of the peculiar binary radio millisecond pulsar 
PSR J1740-5340 and candidate millisecond pulsars (MSPs) in the globular cluster NGC 6397. The 
X-rays from PSR J1740-5340 appear to be non-thermal and exhibit variability at the binary period. 
These properties suggest the presence of a relativistic intrabinary shock formed due to interaction of a 
relativistic rotation-powered pulsar wind and outflow from the unusual "red-straggler/sub-subgiant" 
companion. We find the X-ray source U18 to show similar X-ray and optical properties to those of 
PSR J1740-5340, making it a strong MSP candidate. It exhibits variability on timescales from hours 
to years, also consistent with an intrabinary shock origin of its X-ray emission. The unprecedented 
depth of the X-ray data allows us to conduct a complete census of MSPs in NGC 6397. Based on 
the properties of the present sample of X-ray-detected MSPs in the Galaxy we find that NGC 6397 
probably hosts no more than 6 MSPs. 

Subject headings: globular clusters: general — globular clusters: individual (NGC 6397) — pulsars: 
general — pulsars: individual (PSR J1740-5340) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Globular clusters are well known for th eir abundance 
of ro t ation-powered mi ll isecond pulsars (|Camilo et alJ 
120001 iFreire et all l200l iRansom et all [200l 6 . These 
objects are believed to have been produced by the evo- 
lution of low-mass X-ray binaries (LMXBs) that were 
themselves produced through dynamical interactions of 
binaries and neutron stars. The superb sub-arcsecond 
X-ray imaging capability of the Chandra X-Ray Obser- 
vatory has allowed X-ray studies of MSPs in the dense 
cores of globular clusters for the first time. In addition 
to establishing the X-ray properties of the Galactic pop- 
ulation of MSPs, X-ray studies of cluster MSPs offer con- 
straints on stellar and binary evolution and the internal 
dynamical evolution of globular clusters. To date, X-ray 
counterparts of th ese systems have been detected with 
Chan dra in 47 Tuc jGrindla v et al. 2002: Bogdano v et alJ 
[111 , NGC 6397 ( IGrindla v et al.1120021 and this work), 
M28 (iBecker et al.ll2003l S Bogdanov et al. in prepara- 
tion), M4 (Ba ssa et al 12004), N GC 6752 (iD'Amico et al 
, M71 (|Elsner et al.ll2008f ). and Ter 5 (jHeinke etal 




NGC 6397 is the nearest apparently core-collapsed and 
the second closest globular cluster, with D « 2.4 kpc 
(|Hansen et al.ll200"7UStrickler et al.ll2009l and references 
therein). It hosts one known MSP, PSR J1740-5340 
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(|D'Amico et alJ I2001D , which is bound to an unusual 
"red-straggler" or "sub-subgiant" 7 com panion in a 32.5- 
hour binary orbit (jFerraro et al.l 120031 ). This is at odds 
with the standard MSP formation theory that predicts ei- 
ther a white dwarf or a very- low-mass (~0.03 M^) degen- 
erate companion (see Bhattac harva fc van den Heuvell 
119911 for a review), implying that either this binary has 
j ust emerged from the "recycling" phase in a LMXB 
(jBurderi et al.ll2002l ) or has been involved in a close dy- 
namical binary-binary encounter in which the original 
com panion to the pulsar w as exchanged for the current 
one (|Camilo fc Rasio 20QE and references therein). 

Observations of MSPs, both in the field of the Galaxy 
and in globular clusters, have revealed that their X-ray 
emission can be of thermal and/or non-thermal c har- 
acter (see IZavuril20"07t [Bogdanov & Grindlay HH for 
a review). In many MSPs, the X-rays appear to be 
predominantly due to surface emission from the mag- 
netic polar caps of the neutron star (IBogdanov et al.l 
12001 IZavlinl 12001 IBogdanov fc Grindlavl 12009ft . Non- 
thermal pulsed X-ray emission, on the other hand, can 
arise from particle acceleration processes in the pulsar 
magnetosp here as seen from th e most energetic MSPs 
(see, e.g. iRutledge et al.l [2004). Alternatively, non- 
thermal X-rays can be produced via interaction of the 
rotation-powered pulsar wind with the a mbient interstel- 
lar medium, as seen for PSR B1957+20 (IStappers et all 
2003?) and J2124-3358 (|Hui fc Beckerl l2u06). or material 
from a clo se binary companion a s in PSR J0024-7204W 
in 47 Tuc (|Bogdanov et al.ll2005ft . 

In this paper, we present Chandra deep imaging spec- 
troscopic observations of NGC 6397, with a particular fo- 

7 The terms red straggler and sub-subgiant are used for cluster 
stars that lie in a location below the base of the red giant branch 
that is difficult to understand in the context of standard evolu- 
tion scenarios for single stars or binaries in a dynamically inactive 
environment (Orosz & van Kerkwiik 2003). 
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Chandra ACIS-S 
0.3 - 6 keV 
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Fig. 1. — Coadded Chandra X-ray Observatory 295.1-ks ACIS-S3 image of the core of the globular cluster NGC 6397 in the 0.3—6 keV 
band. The 1.5" circles are centered on the 79 X-ray sources detected within the 2.3' half-mass radius of the cluster (outer dashed circle), 
with the positions of PSR J1740-5340 (U12) and CXOGlb J174041.6-534027 (U18) marked. The inner dashed circle shows the 5.5" core 
radius of the cluster. The grayscale corresponds to number of counts increasing logarithmically from (white) to 3515 (black). North is 
up and east is to the left. 



cus on PSR J1740-5340 and CXOGlb J174041. 6-534027 
(U18), a plausible MSP candidate. We also investigate 
the X-ray source population of NGC 6397 in an attempt 
to identify plausible MSP candidates and constrain the 
MSP content of this cluster. The work is organized as 
follows. In §2 we describe the data reduction and analy- 
sis procedures. In §3 we focus on the properties of PSR 
J 1740-5340, while in §4 we investigate the MSP can- 
didate U18. In §5 we attempt to place interesting con- 
straints on the total number of MSPs in the cluster based 
on the available X-ray and optical data. We present a 
discussion in §6 and conclusions in §7. 

2. DATA REDUCTION AND ANALYSIS 

The Chandra data discussed herein were acquired dur- 
ing two separate observations in Cycle 8, on 2007 June 
22 (ObsID 7461) and 2007 July 16 (ObsID 7460) for 90 
and 160 ks, respectively. In both cases the ACIS S3 chip 



in VFAINT telemetry mode was at the focus. We also 
make use of two Cycle 3 ACIS-S observations of 28.5 ks 
and 27 ks (ObsIDs 2668 and 2669, respectively), and a 
single Cycle 1 ACIS- I (ObsID 79) observation of 49 ks 
(jGrindlav et al.ll2001h . all acquired in FAINT mode. Ta- 
ble 1 summarizes all the observations used in this work. 

The data re-processing, reduction, and analysis were 
performed using CIAO 8 4.0. Starting from the level 
1 data products, we first removed pixel randomization 
from the standard pipeline processing in order to aid in 
source disentanglement in the dense cluster core. In addi- 
tion, for the purposes of faint source detection we applied 
the background cleaning algorithm for the data taken in 
VFAINT observing mode. However, as this procedure 
tends to discard real source counts for relatively bright 

8 Chandra Interactive Analysis of Observations, available at 
htt p : / / cxc . harvard . edu/ciao/ 
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TABLE 1 

Chandra Observations of NGC 6397 



0.5 1 2 5 

Photon energy (keV) 

Fig. 2. — Phase-integrated stacked spectrum of PSR J1740-5340 
from all four Chandra ACIS-S observations, fitted with an absorbed 
powerlaw spectrum. The lower panel shows the best fit residuals. 
See the text and Table 2 for the best fit parameters. 



sources, the background-cleaned data were not used for 
the spectral analyses discussed in §3 and §4 but were 
used for the population study in §5. Before coadding the 
four ACIS-S images, the individual aspect solutions were 
reprojected relative to the longest observation (ObsID 
7460) using the brightest sources in the cluster in order 
to correct for any differences in the absolute astromctry 
between the four exposures. 

X-ray source detection was performed with the 
wavdetect tool in CIAO. In the crowded cluster core, we 
also used the PWDetect 9 script (jDamiani ct al. 1997a1fb1). 
which tends to be more effective at identifying faint 
sources near much brighter sources. The resulting source 
positions were refined using the IDL tool acis_extract. 
The 95% confidence positional uncertainty radius for 
each source was comput ed using the emp irical relation 
given by Equation (5) in lHong et al.l (2005). In total, 79 
sources are detected within the 2.33' half-mass radius of 
NGC 6397 (see Tab l e 3). Compared to the source list 
from IGrindlav et alj (|2001[) , base d on the single ACIS-I 
observation, and IGrindlav! (|2005h . based on the ACIS- 
I and two short ACIS-S observations, we find that all 
sources except U20, U64, U71, U72, U74, U78, and U85 
are detected in the combined deep ACIS-S image. 

For the purposes of optical matching of the X-ray 
sources (see §5) using the set of Hubble Space Telescope 
Advanced Camera for Surveys Wide Field Channel ob- 
servations (GO-10257) in F435W (B), F625W (R), and 
F658N (Ha), we computed a boresight correction for the 
Chandra source coordinates (see H. Cohn, in prepara- 
tion, for details). The resulting shift of the Chandra coor- 
dinates is Aa = -0.19"±0.02" AS = 0.18"±0.02". The 
X-ray position based on this correcti on for PSR J1740- 
5340 agrees with the optical position (jBassa fc Stapperi 
l200l to 0.01" in a and 0.03" in 8. While the Chandra 
ACIS-S mosaic encompasses the entire half-mass radius 
of the cluster, the HST ACS/WFC mosaic provides com- 
plete coverage out to a radius of ~1'.'5 from the cluster 
center, partial coverage in two or more exposures out to 
~2"5, and partial coverage in at least one exposure out 
to 2'.'9. As a result, the sources U5, U16, U77, and U84, 



Telescope/ 


Epoch of 


Observation 


Exposure 


Instrument 


Observation 


ID 


Time (ks) 


Chandra/ACIS-I 


2000 Jul 31 


79 


49.0 


Chandra/ACIS-S 


2002 May 13 


2668 


28.5 


Chandra/ACIS-S 


2002 May 15 


2669 


27.0 


Chandra/ACIS-S 


2007 Jun 22 


7461 


90.0 


Chandra/ACIS-S 


2007 Jul 16 


7460 


160.0 



lie outside the R field making their classification difficult 
due to the abse nce of BR and Ha — R color measure- 
ments. However, iKaluznv et al.l (|2006l ) have shown that 
the optical counterparts of U5 and U77, V30 and V36 
respectively, are variable using ground-b ased photome- 
try a nd thus are probably active binaries (jKaluznv et alJ 
WM). 

Net counts and spectra for each X-ray source were ex- 
tracted using acis_extract from polygonal regions sized 
to contain 90% of the total energy at 1.5 keV. For spec- 
troscopy, the extracted source counts in the 0.3-8 keV 
range were then grouped in energy bins so as to ensure 
at least 15 counts per bin. The X-ray spectral analyses 
of PSR J1740-5340 and U18 were carried out using the 
XSPEC 10 package. For the variability analysis, the pho- 
ton arrival times were first reduced to the solar system 
barycenter, using the CIAO tool axbary. 

3. PSR J1740-5340 

IGrindlav efafl (|2001l . I2002D have found the X-ray 
counterpart of the PSR J1740-5340 system (see Figure 
1) to be a moderately luminous (Lx ~ 10 31 ergs s _1 ) and 
relatively hard X-ray source (r ~ 1.4), compared to most 
MSPs. Indeed, the combined Chandra ACIS-S spectrum 
(Figure 2) is well described by a pure non-thermal model 
with power-law photon index T = 1.73 ± 0.08, hydrogen 
column density of N H = (2.19jlS;||) x 10 21 cm~ 2 , and un- 



9 PWDetect has been developed by scientists at Osservatorio 
Astronomico di Palermo G. S. Vaiana thanks to Italian CNAA 
and MURST (COFIN) grants. 



absorbed X-ray flux F x = (3.2 ± 0.4) x 10~ 14 ergs cm 2 
s _1 (0.3-8 keV) with xl = 1-14 for 42 degrees of free- 
dom. All uncertainties quoted are la. For an assumed 
distance of 2.4 kpc to NGC 6397, the implied X-ray lu- 
minosity in the 0.3-8 keV band is Lx = 2.2 x 10 31 ergs 
s- 1 . 

We also applied a two-temperature thermal (blackbody 
or H atmosphere) model, which provides a good descrip- 
tion of the X-ray spectra of severa l nearby MSPs ( Zavlin 
I2006t Bogdan ov fc Grin dlav 2009). However, this model 
fits the phase-averaged spectrum of PSR J1740-5340 
poorly. The spectrum of PSR J1740-5340 likely con- 
tains a soft thermal component (with T e g ~ 10 6 K and 
Lx of order a fewxlO 30 ergs s _1 ) originating from the 
hot magnetic polar caps of the MSP, as seen in many 
X-ray-detected MS Ps (e.g PSR J0437-4715 and mos t 
MSPs in 47 Tuc, see lZavlin|[2006l : IBogdanov et alj [2006). 
To investigate this possibility, we first introduce a sin- 
gle blackbody component to the spectral model. Fitting 
a composite powerlaw plus blackbody model yields T — 
1.56ljj;& iV H - (2.6±J; 8 8 )xl0 21 cm" 2 , kT = 0.19 +0 04 09 K, 

and Reg = OAbto^l km with xl = 1-1 for 40 degrees of 
freedom. Although an F-test indicates that the addition 

10 Available at http: / /heasarc. nasa.gov/docs/xanadu/xspec/index. html 
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of a thermal component is not statistically warranted by 
the data, the best fit values for T e g and R e g are similar 
to those obtained for the 47 Tuc MSPs so this component 
may in fact be genuine. A composite model consisting of 
a powerlaw plus two thermal components yields good fits 
to the spectrum, although the limited photon statistics 
do not allow meanigful constraints for a three-component 
model. In principle, the presence of surface thermal emis- 
sion could be determined by way of high time-resolution 
X-ray observations, which would reveal any thermal pul- 
sations. Unfortunately, the 3.2-second time resolution of 
the ACIS observations precludes such an investigation 
with the present dataset. 

The X-ray emission from the PSR J1740-5340 system 
could also originate from a thermal plasma within the 
binary, possibly from the active corona of the secondary 
star or the material responsible for the radio eclipses. To 
test this possibility, we used the vmekal thermal plasma 
model in XSPEC, with metal abundance s set to values 
representative of the stars in NGC 6397 (Castilh o et al.l 
l2000t iGratton et al.ll200l . This model also reproduces 
the observed spectral shape rather well, with best fit pa- 
rameters A H = (1.80±S;lf) x 10 21 cm" 2 , kT = 6.2±^ 
keV, and F x = (2.70±g:|§) x 1CT 14 ergs cm 2 s^ 1 (0.3-8 
keV), with \ 2 V — 1-2 for 47 degrees of freedom. As dis- 
cussed in §6.1, although a thermal plasma model is con- 
sistent with the observed X-ray emission from the PSR 
J1740-5340 binary, a predominantly non-thermal origin 
of the X-rays is more likely. 

Based on t h e AC IS-I observation of NGC 6397, 
IGrindlav et al.l (|2002t ) reported evidence for a gradual 
increase (by a factor of 2 in total) in the X-ray count 
rate of PSR J1740-5340. However, due to the lim- 
ited phase coverage and photon statistics, no conclusive 
statements could be made regarding variability. Fig- 
ure 3 shows the X-ray count rate from the pulsar as 
a function of binary phase, based o n the radio timing 
ephemeris from lD'Amico et al.l(|2001[ ). of all Chandra ob- 
servations of PSR J1740-5340. Although the individual 
lightcurves are suggestive of flux modulations, with gen- 
erally lower count rates around the radio eclipse phases 
(0.05-0.45), a Kolmogorov-Smirnov (K-S) test on the 
unbinned lightcurves within each observation reveals no 
statistically significant variability. In addition, the same 
test indicates variability only at a 1.7er confidence level 
over the entire set of ACIS-S observations. On the other 
hand, based on Poisson statistics, we find that the count 
rate minimum that occurs near <fi = 0.25 (see bottom 
panel of Fig. 3) for the lightcurve folded at the binary 
period and grouped in 20-40 bins, deviates by ~3.8<r 
from what is expected from a constant source. A % 2 
test on the same lightcurve indicates a 98.7% probability 
of variability. There is also marginal evidence (at ~95% 
confidence) for spectral variability, with an apparent soft- 
ening of the X-ray emission around <p — 0.25. Due to the 
long binary period, the ACIS-S data cover less than 3 
full binary orbits so it is not clear whether this variabil- 
ity is truly periodic and stable over many orbits. In- 
terestingly, the folded lightcurve exhibits the same char- 
acteristic shape as tha t of PSR J0024-7204W in 47 Tuc 
(Bogdanov et al. 2005), with a minimum around superior 
conjunction, roughly coincident with the radio eclipses. 
The spectral and temporal similarities of the two systems 
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Fig. 3.— Lightcurves of PSR J1740-5340 folded at the pulsar's 
binary period for each Chandra ACIS observation. The bottom 
panel shows the combined folded lightcurve from the four ACIS- 
S observations, with the count rate corrected for the non-uniform 
exposure across the orbit. The phase is defined based on the radio 
pulsar timing convention in which superior conjunction occurs at 
cj> = 0.25. Two orbital cycles are shown for clarity. 



point to an intrabinary shock origin of the X-rays from 
PSR J1740-5340. 

4. U18: A MILLISECOND PULSAR CANDIDATE 

The X-ray source CXOGlb J174041. 6-534027 (here- 
after U18) is positionally coincident with a peculiar "red- 
straggler" star 11 (jGrindlav et al.H2001f) . In addition, it 
is known to have a relatively hard spectrum making it 
a plausible candidate for a MSP given the similarities in 
X-ray and optical properties with PSR J1740-5340. A 

11 We note that the proper motions of the optical counterparts to 
U12 and U18 are consistent with cluster membership (see H. Cohn, 
et al., in preparation). Therefore, their anomalous locations in the 
cluster optical color-magnitude diagram cannot be explained by 
them being foreground or background objects. 
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Photon energy (keV) 

Fig. 4. — Total Chandra ACIS-S spectrum of the X-ray source 
CXOGlb J174041. 6-534027 (U18) fitted with a pure powerlaw 
spectrum. The bottom panel shows the best fit residuals. See 
text and Table 2 for best fit parameters. 



fit of a power-law model to the combined ACIS-S spec- 
trum of U18 gives T = 1.38±g;g|, N H = (1.4±g;^) x 10 21 
cm -2 with xt — 1-0 for 114 degrees of freedom. The 
derived unabsorbed flux of Fx = 9.7 x 10~ 14 ergs cm~ 2 
s- 1 (0.3-8 keV) implies a luminosity of L x = 6.7 x 10 31 
ergs s _1 for D = 2.4 kpc. Figure 4 shows the spectrum 
of U18 in the 0.3-8 keV band and the best fit PL model. 

As with PSR J1740-5340, we also conducted fits using 
a two component (thermal plus non-thermal) model. The 
resulting best fit values are T = l.29t° a % kT = 0.191^ 
keV, R cff = 0.18±°;i 3 , N H = (1.70tg$) x 10 21 cm" 2 , 
with xt = 0.98 for 112 degrees of freedom. The un- 
absorbed flux in the 0.3-8 keV band is 1.07 x 10~ 13 
ergs cm -2 s , which corresponds to an X-ray lumi- 
nosity of L x = 7.4 x 10 31 ergs s _1 for D = 2.4 kpc . 
By analogy with the neares t known MSPs (|Zavlinll2006l : 
iBogdanov fc Grind lav 2009), we also apply a model with 
a powerlaw plus two thermal components. This model is 
consistent with the observed spectrum as well, although 
the relatively high column density along the line of sight, 
which absorbs a large fraction of source photons below 
~0.5 keV, makes it difficult to reliably constrain the 
softer thermal component. 

A thermal plasma model (vmekal with the same as- 
sumed abundances as for PSR J1740-5340) also provides 
a good description of the spectrum of U18. The best 
fit parameters in this case are Nh = (L53lo^ 3 ) x 10 21 
cm- 2 , kT = 16.8±3.I keV, and F x = 9.0 x 10" 14 ergs 
cm~ 2 s" 1 (0.3-8 keV), with xl = 1-06 for 128 degrees of 
freedom. 

It is obvious from Figure 5 that the X-ray flux from 
U18 exhibits variability by as much as a factor of ^4 on 
timcscales spanning from hours to years. Indeed, a K-S 
test on the concatenated set of ACIS-S observations in- 
dicates a 3.5 x 10 -15 probability of the observed photons 
coming from a steady flux distribution. Moreover, the 
photon arrival times from U18 within observations 2668, 
2669, 7460, and 7461 deviate by 3cr, 2.3ct , 3cr, and 4cr 
fr om a constant distrib ution, based on a K-S test. 

iKaluznv et al.l (|2006f ) have reported the likely detec- 
tion of periodicity, possibly due to ellipsoidal variations, 
with P = 1.3 days in the source V31, which they as- 



TABLE 2 

Spectral fits for PSR J1740-5340 and U18 



Parameter 


PSR J1740-5340 
(U12) 


CXOGlb J174041. 7-534027 
(U18) 


Powerlaw 


N n (10 21 cm" 2 ) 

r 

F x a (0.3-8 keV) 
Xl/dof 


2 1Q+ U ' 22 
z ' i9 -0.25 
i 7 o+0.08 
1 - ''-'-0.08 
o 99 +0.37 
°- zz -0.36 
1.14/42 


1 40 +Ui2 
1 - 4U -0.11 
i no+0.05 
1 " 3 °-0.04 
q 72 +0.07 
y - ' z -0.02 
1.00/114 


Powerlaw+Blackbody 


AT H (10 21 cm" 2 ) 

r 

kT (keV) 
R aH h (km) 
F x a (0.3-8 keV) 
Xl/dof 


n ro+l.SU 
z '°°-0.80 

1 56+ 018 
l.OO_ 23 

n IQ+0'09 
u - la -0.04 
15+ 009 
, 47 +0.39 

1.11/40 


-"■•'"-0.07 
-1 9Q +0.05 
1,za -0.12 
n 19+0-03 
u,la -0.03 
0-18_ 14 

10 71+"- 08 
lu ' 1-0.03 

0.98/112 


VMEKAL 


iV H (10 21 cm" 2 ) 
kT (keV) 
F x a (0.3-8 keV) 
x 2 /dof 


1.801";* 
6.2+H 
2 7n+° :29 

z - ' u -0.36 

1.19/47 


, c-q + 0.14 
1 - oo -0.13 

9.04l»; 7 I 
1.06/128 


a Unabsorbed flux in 
keV band. 

b Effective blackbody 
kpc to NGC 6397. 


units of 10 14 ergs cm 2 s 1 in the 0.3—8 
emission radius assuming a distance of 2.4 



sociate with U18. It should be noted, however, that 
the position of V31 lies at ~2.1<t from U18 after bore- 
sight correction. Moreover, V31 does not correspond to 
the red-straggl e r opti cal counterpart of U18 reported by 
IGrindlav et al.l (|2001l ), which is only ~0.4er away. We 
have folded the set of four Chandra ACIS-S observations 
at the period of V31 but find no statistically significant 
variability. This implies that V31 is most likely not the 
actual optical counterpart of U18. 

5. AN X-RAY CENSUS OF MILLISECOND 
PULSARS IN NGC 6397 

As noted previo usly, several nearby MSPs in the 
field of the Galaxy (jZavlinl 120061: IBogdanov fc Grindlavl 
12009ft and most MSP s in the globular cluster 47 Tuc 
(Bog danov et al.ll2006l ) appear to be predominantly soft, 
thermal X-ray s ources due to their hot polar caps. As 
demonstrated by IBogdanov et akl (|2008D (see in particu- 
lar their Figure 1), the (nearly) antipodal geometry of the 
polar caps and the effects of light bending guarantee that 
this surface emission is observed for any combination of 
magnetic inclinations and viewing angles. In addition, 
the difference in observed luminosity between the most 
and least favorable geometric orientations (always face- 
on versus always edge-on hot spots) is only a factor of ~3 
as a consequence of gravitational bending of light near 
the neutron star surface. Thus, thermal MSPs should 
be observable in X-rays even if they cannot be seen in 
the radio due to unfavorable beaming and fall within a 
relatively narrow range of X-ray luminosities. The limit- 
ing sensitivity of the available ACIS-S data is ~10 29 ergs 
s , meaning that typical thermal MSPs (with luminosi- 
ties of a few xlO 30 ergs s _1 ) would be detected in the 
combined Chandra ACIS-S observations if they were lo- 
cated in NGC 6397, even in the unlikely event that their 
two hot spots are always seen edge-on. Binary MSPs 
that eclipse in the radio tend to have X-ray luminosities 
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Fig. 5.— X-ray lightcurves of CXOGlb J (U18) in the 0.3-8 
keV band from all Chandra ACIS observations binned in 1-hour 
intervals. Time zero corresponds to the start of each observation. 



~10 3 1 ergs s" 1 (IStappers et all l2003t iBogdanov eFall 
[20051 [200a lElsner et all 120081 ). Moreover, the intrabi- 
nary shock in such systems should presumably be visible 
in X-rays regardless of the pulsar geometry and the ori- 
entation of the binary. Therefore, based on our current 
understanding of the X-ray properties of these sources, 
all MSPs in NGC 6397 should be detected with high 
significance in the available data. 

We have compiled a set of "template" MSP X-ray col- 
ors and luminosities using the cluster and field MSPs that 
have been detected in X-rays t o date. These includ e the 
11 MSPs listed in Table 1 of iPavlov et~aT] (|2007l ). the 
15 MSPs in t he globular clusters 4 7 Tuc with uncon- 
fused spectra (IBogdanov et all [2006). one MSP in M71 
(lElsner et all 120081 ). two MSPs in Ter 5 (jHeinke et all 
|2006J), and the recentl y discovered PSR J10 23+0038 in 
the field of the Galaxy (| Archibald et al.ll2009f ). Using the 
measured spectral properties of these MSP we calculate 
their X-ray luminosities and colors as they would appear 
in NGC 6397 using N H = 2 x 10 21 cuT 2 and D = 2.4 




-2 2 -2 2 

Xcolor=2.5 log([0.5-1.5 keV]/[1.5-6 keV]) 

Fig. 6. — (a) X-ray luminosity versus color diagram for NGC 
6397. The black points show all 79 sources within the NGC 6397 
half- mass radius (Table 3), for an assumed powerlaw spectrum with 
r = 2.5. (6) Same as (a) but with all sources with probable optical 
identifications removed. The red circles in both panels correspond 
to simulated MSPs based on the current Galactic sample of MSPs 
detected in X-rays (see text for details). Sources with no optical 
counterparts arc marked with black squares, those with ambiguous 
or unknown optical counterparts with solid blue triangles, the main 
sequence star candidates with black crosses, and the main-sequence 
turnoff star candidate with the open magenta circle. The dotted 
lines extending below the error bars for the MSPs illustrate the 
effect of the uncertainty in the viewing geometry of the X-ray- 
emmiting polar caps. 



kpc and compare them to the cluster X-ray source pop- 
ulation (see Figure 6). We define the color as Xcolor = 
2.51og(counts[0.3 - 1.5 keV]) /counts [1.5 - 6 keV]) and 
consider the absorbed luminosity in the 0.5-6 keV band. 
To estimate the uncertainties in the derived MSP lumi- 
nosities, we have considered a nominal 30% error in the 
pulsar distances plus the uncertainty due to the range of 
possible distances to NGC 6397. Furthermore, we have 
taken into account the lack of information regarding the 
viewing angles of the X-ray-emitting polar caps of any 
thermal MSPs in NGC 6397 (represented by the dotted 
lines in Figure 6). High-quality X-ray spectra of nearby 
thermal MSPs show at least two thermal components 
and a much fainter hard tail above ^2 keV. For crude 
photon statistics, the X-ray emission of these MSPs re- 
sembles a powerlaw with photon index r ~ 2 — 3 in the 
0.5-6 keV range. Thus, to estimate the X-ray luminosi- 
ties of the sources in NGC 6397 we consider T = 2.5. 
To ascertain the effect of our choice of spectral model 
we also used a pure blackbody with kT — 0.2 keV and 
varied AT H = (1 - 3) X 10 21 cm" 1 . We find that the as- 
sumptions regarding the spectral shape do not affect the 
conclusions of our analysis. 

Using the computed set of MSPs, we attempt to iden- 
tify MSPs in NGC 639 7 without the be nefit of ra- 
dio timing searches (see IHeinke et al.l 12001 for a sim- 
ilar study for 47 Tuc). Important constraints can be 
gained from optical counterpart identifications using 
HST ACS/WFC observations of the core of NGC 6397 
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(H. Cohn et al. in preparation). Of the 79 X-ray 
sources detected within the half-mass radius of NGC 
6397, 64 sources have probable associations in the op- 
tical (see Table 3) predominantly with active binaries 
(ABs) and cataclysmic variables (CVs). Additionally, 
U24 is classified as a qLM XB based on its X-ray proper- 
ties (|Grindlav et al.ll200l . The remaining sources either 
have ambiguous/unknown or no optical counterparts. 

Solitary MSPs should have virtually no optical coun- 
terparts 12 . Therefore, any X-ray source without an 
optical counterpart is a viable MSP candidate. Bi- 
nary MSPs both in clusters and the field are com- 
monly bound to low-mass He white dwarfs or very-low- 
mass (brown-dwarf-like) companions (~0.03 M Q ) such as 
in the canonical " bla ck- widow" system PSR B1957+20 
( Fructiter et all fl988h . For NGC 6397, IStrickler et all 
(2009) have found that none of the 41 He white dwarfs 
identified optically have X-ray counterparts, implying 
that they are not partnered with MSPs. The optical 
properties of black widow systems are not well estab- 
lished as only two obj ects have been studie d in any de- 
tail, PSRs B1957+20 dRevnolds et al.l l200l and J2051- 
0827 (Sta ppers et alJ 120011 ). Nonetheless, these systems 
tend to have Ly > 10 31 ergs s" 1 (IStappers et "ail 120031: 
iBogdanov et"aTll2006t lElsner et all 120081 ) so thelack of 
unidentified sources with comparable X-ray luminosities 
implies the absence of black widow systems in NGC 6397. 
Moreover, this indicates that aside from J1740-5340 and 
U18, there are no additional peculiar interacting bina- 
ries in this cluster either. Only a single source without 
an optical counterpart (U113 in Table 3) has X-ray colors 
and luminosity (Xcolor = 1.05 and logLx = 30.5) con- 
sistent with those of a typical MSP (Figure 6). The sole 
source with an ambiguous/unknown optical counterpart 
that has soft colors, U16, as well as U41 and U112, which 
appears to hav probable main-sequence optical counter- 
parts, are also possible MSP candidates. Based on these 
findings, we obtain a limit of <6 MSPs in NGC 6397. 

6. DISCUSSION 
6.1. PSR J1740-5340 

Optical observations of the PSR J 1740-5340 system 
(|Ferraro et al.ll2003t ISabbi et al.ll2003[ ) have revealed the 
presence of ellipsoidal variations as well as a stream of 
gas protruding from the 1.6 Rq secondary star through 
the LI point implying that the companion is Roche-lobe 
overflowing. However, the fact that PSR J1740-5340 is 
observed as a radio pulsar implies that this gas stream 
never reaches the underlying NS and is instead swept 
back and probably expelled from the binary system by 
the relativistic particle wind from the MSP. The presence 
of this swept-back gas stream i s suggested by th e very 
unusual Ha emission line profile (jSabbi et al.ll2003fl . This 
gas is likely the cause for the irregular radio eclipses the 
MSP exhibits over a wide range of orbit al phases. 

As seen in "black- widow" systems (Stappers et al.l 
120031: IBogdanov et al.ll2"00l lElsner et al.ll2008D and "ex- 
changed" MSP systems (|Bogdanov et al.l 120051 ) the in- 

12 The current best limits on optical emission from isolated 
MSPs come from observations of the nearby PSRs J0030+0451 
(D « 300 pc) and J2124-3358 (D fa 270 pc), for which no coun- 
terparts are found down to B, V, and R magni tudes of 27 — 28 
(Ko ptsevich et al.l|2u~0 3; Mignani & Becker 2O03D- 



teraction between the pulsar wind and material from the 
companion star should result in a shock front where the 
ram pressure of the wind balances the pressure of the in- 
falling gas leading to X-ray emis sion due to particle ac- 
celeration (jArons fe Tavanilll993f ). In this scenario, syn- 
chrotron emission is expected to be the primary energy 
loss mechanism in the resulting shock wave, given that 
it occurs in a relatively strong magnetic field. In the 
PSR J1740-5340 binary, for a shock situated approxi- 
mately at the LI point and an isotropic pulsar wind, the 
pulsar magnetic field at this distance (sa4 x 10 11 cm) 
immediately upstream from the shock is B\ w 2 — 3 
G or 5i w 0.1 G , implying a po st-shock field of 
B 2 = 3£?i ~ 6 - 9 G or 5 2 ~ 0.3 G (|Arons fc Tavanil 
1993). The two values correspond to the two possible 
cases of a magnetically dominated shock (cr ^> 1, where 
a is the magnetic to kinetic energy flux ratio) or a ki- 
netic energy dominated shock (a = 0.003 as in the Crab 
pulsar). This implies that the X-ray emission from PSR 
J1740-5340 source is most likely non-thermal in origin 
and not from a thermal plasma. 

An intrabinary shock origin of the X-rays from 
J1740-5340 is also suggested by the fact that the de- 
rived spectral parameters are quite similar to those 
of the canonical "bla ck-widow" pulsar, PSR B1957+20 
(|Stappers et al.ll2003f) . the peculiar MSP -mam-sequence 
binary in 47 Tuc, PSR J0 024-7204W (jEdmonds et all 
I2002L Boadano v et al.l 12005). and the recently identified 
radio MSP in the field of the Galaxy, PSR J1023+0038 
( Archibal d et alJ 120091 and references therein), which 
also appears to be bound to a non-degenerate ("main- 
sequence-like") star. Moreover, the shape of the X-ray 
modulations at the binary period are quite similar to tha t 
of J0024-7204W (see Figure 1 in IBogdanov et alJl2005Q . 
This variability can be attributed to an occultation of the 
shocked material by the secondary star. This scenario is 
supported by the compelling (albeit marginal) evidence 
for a softening of the spectrum at <j> ~ 0. 25. The con- 
strain t on the system inclination (47° — 56° iFerraro et al.l 
2003:) suggests that the line of sight to the pulsar is never 
obstructed by the companion. Therefore, when the shock 
is occulted at <f> sa 0.25 (either partially or totally), the 
soft thermal radiation from the MSP surface provides a 
larger contribution to the total X-ray emission, resulting 
in a softening of the observed spectrum. In principle, 
detailed orbital-phase-resolved X-ray observations would 
permit disentanglement of the thermal polar cap emis- 
sion and the non-thermal shock emission. 

Note that the moderate spin-down luminosity 13 of 
E sa 3.3 x 10 34 ergs s" 1 (|Bassa fc Stappers! l200l does 
not favor a magnetospheric origin of the observed non- 
thermal X-rays considering that most other MSPs with 
comparable values of E exhibit much softer, predomi- 
nantly thermal spectra (see, e.g.. IBogdanov et all 120061 : 
lZavliijl2004 IBogdanov fc Grindlavll2009f ). Furthermore, 
for pulsed non-thermal emission, it is difficult to explain 
the observed variations in the total flux and the softening 
of the spectrum at <f> s=a 0.25. 

6.2. U18: A Hidden MSP? 

13 As discussed by)D'Amico et al. (2001), the effect of cluster 
acceleration on the measured pulsar spindown rate is neglegible. 
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The X-ray and optical properties of U18 are consistent 
with those of a binary containing a rotation-powered pul- 
sar wind interacting with material from the secondary 
star. It is interesting to note the significantly larger lu- 
minosity compared to J 1740-5340 and other interacting 
MSP systems. This could indicate either a much larger 
shock region, an enhanced density of radiating particles, 
and/or a more energetic pulsar wind. If this binary does 
indeed harbor a MSP, it may be difficult to detect at ra- 
dio frequencies due to the large quantities of gas present 
within and around the binary. This gas may render the 
pulsar perpetually eclipsed at radio frequencies. Thus, 
U18 may belon g to the class of so-called "hidden" MSPs 
(|Tavani 11190 lh making a confirmation of its true nature 
quite difficult. The detection of pulsed X-ray emission 
from the pulsar itself may also be difficult due to the 
dominant X-ray flux from the intrabinary shock and the 
a priori unknown pulsar spin period. 

We note that, at present, we cannot strictly rule out 
the possibility that UI8 is a cataclysmic variable, such as 
an intermediate polar (IP), instead of an MSP. These ob- 
jects also generally exhibit hard X-ray spectra and X-ray 
luminosities compar able to that of U18 . In addition, the 
IP AK09 in 47 Tuc (jHeinke et al.ll2005ft appears to have 
a "red straggler" or "sub-subgiant" companion as well. 
In order to unveil the true nature of the UI8 system, fur- 
ther optical observations, optimized for variability study, 
are needed. The cadence of our HST ACS/WFC optical 
dataset (10 single-orbit observations spaced ~1 month 
apart, see H. Cohn et al. 2009, in preparation) is not 
suitable for a reliable determination of the binary period 
of this system. 

6.3. The Millisecond Pulsar Population of NGC 6397 

The analysis in §5 suggests that NGC 6397 most 
likely harbors 1 — 6 MSPs. It is of interest to com- 
pare these limits to the expected number of MSPs in 
this cluster, scaling by the stellar encounter rate (r, 
iVerbunt fc Hutl 119871) from oth er clusters with detected 
MSPs 14 . iPoolev et all (|2003h integrate collision rates 
over King models out to the half-mass radius, finding 
a collision rat e for NGC 6397 1/74 that of 47 Tuc. O ther 
analyses (e.g. IHeinke et alj|2003t iBassa et alj 12008ft . us- 
ing a simpler formalism (r = i"c> where po is the cen- 
tral density and r c is the core radius of the cluster) , find 
collision rates for NGC 6397 roughly 1/10 that of 47 Tuc. 
As 47 Tuc has 23 known MSPs, and has been sugg ested 
to have a total of -25 MSPs (jHeinke et all 12005ft . the 
two approaches suggest either < 1 or ~2 MSPs should be 
present in NGC 6397. Comparisons with the number of 
MSP s in other clusters, such as Terzan 5 (Ra nsom et al.l 
2005) give similar results. However, T comparisons seem 
to un derpredict the num ber of X-ray sources in NGC 
6397 (jPoolev et al.l [2003ft . indicating differences in the 
X-ray binary production processes in cor e-collapsed clus - 
ters such as NGC 6397. For instance, iFregeaul (|2008ft 
has suggested that continuing globular cluster collapse 
means that the values of T of King-model clusters dur- 
ing the epoch of X-ray binary formation were a factor of 

14 We note that the numbers of (quiescent) LMXBs in globu- 
lar clusters, likely progenitors of MSPs, appear to scale with V 



3 smaller. In any case, the predicted numbers of MSPs 
in NGC 6397 from any of the above comparisons are gen- 
erally in agreement with our estimate of 1 — 6 MSPs in 
NGC 6397. 

7. CONCLUSION 

We have presented Chandra X-ray Observatory ACIS- 
S observations of the nearby core-collapsed globular clus- 
ter NGC 6397. The depth of the available data has al- 
lowed interesting insight into the MSP population of this 
cluster. For instance, it has provided a more complete 
multi-wavelength picture of the PSR J1740-5340 binary. 
This peculiar system appears to exhibit predominantly 
non-thermal emission modulated at the orbital phase. 
Given the pulsar energetics, the variable non-thermal X- 
ray emission is indicative of interaction between the rel- 
ativistic pulsar wind and material from the companion, 
with evidence for an occultation of the resulting shock 
by the secondary star. Ascertaining the detailed geom- 
etry of the intrabinary shock would require much more 
sensitive orbital-phase resolved spectroscopic X-ray and 
optical observations of this system than presently avail- 
able. 

Based on its X-ray and optical properties, U18 is a 
strong candidate for a MSP, though one that may be 
perpetu a lly ec lipsed at radio frequencies. As noted by 
iTavani I (|T991) such MSPs may be relatively common 
both in globular clusters and in the field of Galaxy More- 
over, these sources may be mis-classified as LMXBs or 
CVs based o n their X-ray to opt i cal flux ratio alone (se e 
in particular [Homer et al.l 120061 lArchibald et al.l [2009). 
Further detailed multi-wavelength observations of PSR 
J1740-5340, U18, and similar systems may reveal more 
information concerning the physics of the shock, which, 
may, in principle offer insight into the properties of MSP 
winds, collisionless relativistic shocks, and particle ac- 
celeration mechanisms. Moreover, this and similar MSP 
binaries provide important clues about the behavior of 
accreting neutron stars transitioning to rotation power, 
which are believed to also contai n an active radio pul- 
sar and an overflowing co mpanion (jCampana et al.lf2004h 
IHeinke et al.l [20071 [2009ft . 

The deep X-ray observations of NGC 6397 supple- 
mented by optical observations, also provide insight into 
the MSP population of this cluster. Specifically, if the 
X-ray luminosity function of the NGC 6397 MSP popula- 
tion resembles that of the Galactic sample of MSPs, then 
this cluster contains no more than —6 MSPs. This range 
of values is generally consistent with those inferred based 
on scaling of stellar encounter rates in globular clusters. 
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a Source names based on the X-ray source list from lGrindlav et al.l pOOll l. 

b J2000.0 acis_extract derived positions of each source bores i ghted to the frame of the HST ACS dataset. 
c The 95% positional uncertainty radius based on I Hong et alj (2005), in units of arcseconds. 
d Background-subtracted source counts in the 0.3-1.5 and 1.5-6 keV bands. 

e Source flux in the 0.5—6 keV band in units of 10 — 16 ergs cm -2 s _1 assuming a power-law model with spectral photon index T = 2.5. 
* Probable source identifications based on optical counterpart properties (from H. Cohn et al., in preparation). AB is an active binary, CV 
a cataclysmic variable, GLX a background galaxy, qLMXB a quiescent low mass X-ray binary, MS a main sequence star, MSTO a main 
sequence turn-off star, and MSP a millisecond pulsar. Sources marked with an additional "?" are probable but less certain classifications, 
those with only "?" have ambiguous or unknown, and those labeled "none" have no optical counterparts. 



